Mesoscopic electron and phonon transport through a curved wire 
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There is great interest in the development of novel nanomachines that use charge, spin, or energy 
transport, to enable new sensors with unprecedented measurement capabilities. Electrical and ther- 
mal transport in these mesoscopic systems typically involves wave propagation through a nanoscale 
geometry such as a quantum wire. In this paper we present a general theoretical technique to de- 
scribe wave propagation through a curved wire of uniform cross-section and lying in a plane, but 
of otherwise arbitrary shape. The method consists of (i) introducing a local orthogonal coordinate 
system, the arclength and two locally perpendicular coordinate axes, dictated by the shape of the 
wire; (ii) rewriting the wave equation of interest in this system; (Hi) identifying an effective scatter- 
ing potential caused by the local curvature; and (iv), solving the associated Lippmann-Schwinger 
equation for the scattering matrix. We carry out this procedure in detail for the scalar Helmholtz 
equation with both hard-wall and stress-free boundary conditions, appropriate for the mesoscopic 
transport of electrons and (scalar) phonons. A novel aspect of the phonon case is that the reflec- 
tion probability always vanishes in the long-wavelength limit, allowing a simple perturbative (Born 
approximation) treatment at low energies. Our results show that, in contrast to charge transport, 
curvature only barely suppresses thermal transport, even for sharply bent wires, at least within the 
two-dimensional scalar phonon model considered. Applications to experiments are also discussed. 

PACS numbers: 85.85.+j, 03.65.Nk, 63.22.+m, 73.23.-b 



I. INTRODUCTION 

A new class of nanomachines is attempting to mea- 
sure extremely minute amounts of energy, of the or- 
der a few neV, and to use such calorimeters to probe 
fundamental properties of thermal conduction in the 
nanoscale regime iSii 4 ^ Like the related case of electri- 
cal conduction;^ low-temperature thermal conduction 
in nanostructures is entirely different than in macro- 
scopic materials because the phonons are in the meso- 
scopic regime, where they scatter elastically but not in- 
elastically. Because inelastic scattering is required to 
establish thermodynamic equilibrium, there is a break- 
down of Fourier's law and the heat equation, which 
assume a local thermodynamic equilibrium character- 
ized by a spatially varying temperature profile. These 
nanodevices have inspired considerable theoretical work 
on thermal transport by phonons in the mesoscopic 
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In this paper we introduce a general method to calcu- 
late the scattering matrix for waves propagating through 
a curved wire or waveguide. The wire is assumed to be of 
uniform cross-section and lying in a plane, but the curved 
segment may have any smooth curvature profile^ such 
as that shown schematically in Fig. ^ The ends of the 
wire (the "leads") are also assumed to be straight. For 
defmitcness we consider two-dimensional waves described 
by the scalar Helmholtz equation 
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FIG. 1: An example of the type of two-dimensional curved 
wire geometry considered in this paper. The shape of the wire 
is used to dehne local orthogonal coordinates X and Y , with 
X the arclength along one of the edges. The width of the wire 
is b. 



of energy e and bulk sound velocity v. Electron spin is 
of no importance here and is neglected. The boundary 
conditions at the edges of the wire are 



$(r) 
n • V$(r) 



0, 
0, 



(for electrons) 
(for phonons) 



(2) 



[V 2 +a]$(r) = 0, r=(x,y) 



(1) 



which is appropriate for electrons or scalar phonons in 
flat wires with rectangular cross-sectioni^i Here a(e) = 
2me/h 2 in the case of electrons of energy e and mass m, 
whereas a(e) = e 2 /h 2 v 2 in the case of scalar phonons 



where n(r) is a local outward-pointing normal. Here we 
have assumed conventional hard-wall boundary condi- 
tions for the electronic states, but stress-free conditions 
for the elastic waves because in this case the wires are 
usually freely suspended. 

Although the wave equation in Eq. is certainly sim- 
ple, the scattering problem described here is complicated 
because the boundary conditions are applied along 
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the curves defining the two edges of the wire. Our ap- 
proach involves rewriting Eqs. (0) and J2} in terms of 
new curvilinear coordinates X and Y, dictated by the 
shape of the wire, such that the wave equation becomes 
more complicated (the wire's curvature produces an ef- 
fective potential), but the boundary conditions become 
trivial. We choose X to be the arclength along one edge 
of the wire, and Y is locally perpendicular. This trans- 
formation allows us to apply the standard techniques of 
scattering theory, including solution of the Lippmann- 
Schwingcr equation, in the XY frame. 

A particularly novel aspect of the phonon transport 
problem is that the reflection probability always van- 
ishes in the long- wavelength limit, permitting an analytic 
(second-order Born approximation) treatment at low en- 
ergies. The energy- dependent transmission probability is 
then expressed simple functional of the curvature 

profile k(X), making possible a straightforward analysis 
of a variety of wire shapes. The fact that long-wavelength 
phonons have perfect transmission is a consequence of the 
rigid-body nature of the underlying system: An clastic 
wave with infinite wavelength is just an adiabatic rigid 
translation of the wire, which must transmit energy per- 
fectly. 

There has been considerable attention given to meso- 
scopic electron transport through curved wires and 
waveguides^ but none to thermal transport. Electron 
transmission probabilities in curved wires are usually ob- 
tained by mode-matching, a method restricted to piece- 
wise separable geometries (wires composed of straight 
segments, circles, and other shapes where the wave equa- 
tion is separable). A related problem that has been stud- 
ied extensively is the formation of electronic bound states 
and resonances in curved wires, where the mapping to lo- 
cal curvilinear coordinates is also often used^^^ Sur- 
prisingly, we are not aware of any work using moving 
frames and then directly solving the resulting Lippmann- 
Schwingcr equation in that basis^& Nor are we aware of 
the use of this method in the extensive microwave engi- 
neering literature^SSiffi where the (more generally appli- 
cable but purely numerical) finite-element method is the 
technique of choice. 

In the next section we carry out the above analysis for 
the two-dimensional Helmholtz equation. In Sec. IIIII we 
consider electron transport through a circular right-angle 
bend, recovering results obtained by Sols and Macucci^i 
and by Lin and Jaffe2£ using mode-matching methods. 
Our main results are given in Sec. II VI where we address 
thermal transport through curved wires. Sec. IVl contains 
a discussion of our conclusions and the experimental im- 
plications of this work. 



II. APPLICATION TO SCALAR WAVE 
EQUATION 



A. Curvilinear coordinate system 

First we use the shape of the wire to define a curvilin- 
ear coordinate system, the arclength X along one edge 
and a locally perpendicular coordinate Y. Which edge 
one chooses is of course arbitrary. The direction of in- 
creasing Y will be chosen so that X, Y, and z form a 
right-handed coordinate system. Both edges are assumed 
to be a smooth plane curves^ 

The unit tangent vector 



e.x 



dr 

dX 



dx 
dX 



dy 
dX 



and normal 



ey 



dy 
' dX 



dx 
dX 



(3) 



(4) 



define local orthonormal basis vectors for the XY frame. 
The sign in Eq. |@J is chosen so that ex x ey = e 2 . We 
then use the Frenet-Serret equation 



de 



x 



dX 



k(X) ey 



(5) 



to define a signed curvature k(X) of the Y = edge. 
According to this definition, k(X) is positive when Y 
increases toward the center of curvature. We will also 
make use of the metric tensor 



(1-kY) 2 
1 



(6) 



in the XY system. 



B. Helmholtz equation in curvilinear coordinates 

Next we make a coordinate transformation from r to 
R = (X, Y), and rewrite the wave equation in terms of 
these coordinates. A convenient way to do this is to use 
the identity 

V 2 = (det 5 )-^(det,)^-^. (7) 
The Helmholtz equation (JTJ then becomes 

[a| + a^-y + a]$(R) = o, (8) 



d Y from 



where we have separated the combination d\ 
the many terms that appear on the right-hand-side of 
Eq. 10), and combined the remaining ones into an effec- 
tive potential V. The potential V is itself a differential 
operator; an explicit expression will be given below. The 
boundary conditions J5J are now 



We now explain our method in detail and apply it to 
the scalar scattering problem stated in Eqs. (0) and J2J. 



<5(R) = 0, (for electrons) 
<9y$(R) = 0, (for phonons) (9) 



200 



FIG. 2: Scattering problem in the XY frame. Here the 
wire appears straight, but the curvature induces an effective 
potential V that causes scattering. 



on the edges Y = and Y = b, with b the width of the 
wire. A scattering state $(R) becomes fully determined 
once we specify its behavior as X — > ±oo. 

In the XY frame the wire appears straight, as illus- 
trated in Fig. [3 The scattering potential vanishes in the 
leads because the wire is straight there. We are now able 
to use conventional scattering theory. 



C. Unperturbed scattering states 

The unperturbed (V = 0) scattering states for both 
spinlcss electrons and scalar phonons arc labeled by three 
quantum numbers cr, n, and e, and can be written in the 
form 



0crne(R) 



Xn(Y), 



(10) 



(11) 



where 

kne = V a{e) - (nir/b) 2 
is the wave number along the wire, 

Xn( Y)= i VVbMmr Y/b) 

1 \J (2—S n o)/b cos(mrY/b) for phonons 

(12) 

is a trigonometric function satisfying the transverse 
boundary conditions of Eq. © , and 



for electrons 



1 



27T 



dk ne 



de 



1/2 



(13) 



is a real normalization constant, cr is a chirality index, 
defined by 



+ 1 if moving in +X direction 
— 1 if moving in —X direction 



(14) 



and n is an integer-valued branch index. The transverse 
eigenfunctions Xn{Y) are normalized according to 



dY X n(Y)xn>(Y) = 5 n 



(15) 




FIG. 3: Dispersion relation for unperturbed scattering states, 
for both spinless electrons and scalar phonons. ot is equal to 
2me/h 2 in the case of electrons, or e 2 /h 2 v 2 in the case of 
phonons. a is 1 on the right half of the figure, and —1 on the 
left, b is the wire width. 



The dispersion relation given by Eq. I|ll|) is shown in 
Fig. El 

The allowed values of the quantum numbers cr, n, and 
e are as follows: The allowed energies form a continuum, 
from e m ; n to oo. Here 



t tz2 2 

_ I f° r electrons 







for phonons 



(16) 



For each value of energy, the branch index takes the val- 
ues in a set S of integers defined by 

g _ \ 1, 2, . . . , N for electrons 
1 0, 1, 2, . . . , N for phonons 

where 



7V( £ )^e[a(e)-^] - 1, 



(18) 



n=0 



with Q(x) the unit step function. For electrons, N is the 
number of propagating channels below energy e, whereas 
for phonons the number of propagating channels is N+l. 
Finally, for each allowed value of e and n, a takes on the 
values ±1. 

The free scattering states satisfy the orthonormality 
and completeness conditions 

J d 2 R ^„ e (R) ^We'(R) - <W <W S(e - e') (19) 

and 

de E E ^neW <WR') = S(R - R')> (20) 

-min a =±lneS 

where n takes the values given in Eq. I|17[l. 
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D. Effective potential 

The terms d x + d\ have been separated out from 
the right-hand-side of Eq. Q so that the eigenfunctions 
<£>(R) reduce to that of a straight wire when V = 0. Ac- 
cordingly, the effective potential is given by 



V 



k'Y 



pP- 

(1 - kYY x (1 - kY)3 ~~ ' 1 - kY 



d 



'x- 



Thcre will be no singularities in V as long as 
— oo <«£><!. 



dy. (21) 



(22) 



The condition l|22|l guarantees that both the Y = and 
Y = b edges of the wire are smooth. 

In applications where the radius of curvature is 
much larger than b, Eq. (|21|l can be simplified. To leading 
order in Kb the effective potential reduces to 



V = -2kY3 2 x - KY(l + 3KY)d x + Kd Y - 



E. Lippmann-Schwinger equation 



(23) 



The scattering problem in the XY frame can be solved 
by standard methods. The Lippmann-Schwinger equa- 
tion for an eigenfunction ^(R) with (electron or phonon) 
energy e is 

$(R) = & n (R) + / d 2 R> G (R, R', c) V *(R'), (24) 

where ^in(R) is a free scattering state coming in from 
the left, and where Go(R, R',e) is the Green's function 
for the unperturbed Hclmholtz equation, satisfying 



[d 2 x + dl + a(e)] G (R, R', e) = 5(R - R) 



(25) 



along with the boundary condition that Go(R, R', e) van- 
ishes as \X — X'\ — > oo. The Lippmann-Schwinger equa- 
tion gives the solution of Eq. JSJ subject to the condi- 
tion that $(R) reduces to the incident state </>; n (R) when 
X -> -oo. 

The unperturbed Green's function for both electrons 
and phonons is 



Gq(R, R', e) 



*_ Xn(Y)Xn(Y ) iknAX -X'\ 

2 ^ k 



(26) 



We note that the summation in Eq. (|26|l is not restricted 
to the values given in Eq. (|17|l . In particular, off-shell 
values of k ne are included. Furthermore, in the electron 
case the n = term in the summation vanishes, because 
the transverse eigenfunction XoO^) vanishes. 

We will also need to write Eq. (|24[l in the alternative 
form 

$(R) = in (R) + J d 2 R' G(R, R, e) U0 in (R), (27) 



where G(R, R', e) is the Green's function of the perturbed 
Helmholtz equation, satisfying 

G(R,R',e) = G (R,R',e) 

+ J d 2 i?"G (R,R",e)FG(R",R',e). (28) 

F. Transmission probability 

The transmission coefficient matrix t nn t(e) gives the 
probability amplitude for a right-moving electron or 
phonon to forward scatter from branch n to branch n' 
at energy e. We define t nn i (e) to be zero if one or both 
branches have minima above e. 

We can obtain a formal expression for t nn i (e) by writ- 
ing the X — > oo limit of the unperturbed Green's function 
as 



Go(R, R', e) 



ly. 

2 ^ 



s (R)^„ e (R') 



(29) 



where the subscripts R denote right-moving (a = +1) 
waves. The summation in Eq. (|2*9")> is now restricted to 
n < N because the higher lying contributions are expo- 
nentially small in the X — > oo limit. Then from Eqs. i|24|) 
and 12911 we obtain 



lim *(R) = <W(R) t £ £ 

X— >00 2 Z ' 



9 h r 2 



(30) 

where the incoming right-moving wave is assumed to be 
in channel n\. Therefore we conclude that 



lim $ = 



'Rme 



and 



where 



lim $ = y t ni n(e) 4>Rne, 



tnn'\£) — Own' o 7 2 



(31) 



(32) 



(33) 



is the amplitude for an incident wave (R, n, e) to forward- 
scatter to (R, n',e). tnn'(e) is called the transmission 
matrix. 

We emphasize that for the case of electron transport, 
t(e) is an N x N matrix, where N varies with energy 
as indicated in Eq. (|18|) . For phonons, t(e) is 2V+1- 
dimensional. 

The expectation value in Eq. I|33|) involves the exact 
scattering state $ with boundary condition correspond- 
ing to an incident state (R, n, to). Using Eq. H27J) we write 
this as 

(<fa.n>z\V\§) = {4>K. n >e\V\<fa.ne) 

+ J <^ n , e (R) V R G(R, R, e) V R . ?W(R'), (34) 



FIG. 4: Section of quantum wire with circular right-angle 
bend. 
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e/A 



10 



12 



where V R and V R i act on R and R', respectively. The 
numerical method we use consists of expressing Eq. <|28[) 
in a basis of unperturbed scattering states, solving this 
equation by matrix inversion, and using Eq. I(34|) to ob- 
tain the transmission matrix in Eq. i|33|) . 

Finally, we define the energy-dependent transmission 
probabilities T e i and T p h for electrons and phonons that 
determine the electrical and thermal currents. For elec- 
trons, the relevant quantity is the ratio of transmitted to 
incident charge current, given by 



Tei(e) 



N 



= Trt+i 



(35) 



For thermal transport by phonons, the relevant quantity 
is the fraction of transmitted energy current given by 



Tph(e) 



N 

E 

71,77,' —0 



(36) 



where v„(e) is the phonon group velocity in the straight 
wire, which can be written as Hv 2 k ne /e. Because the bulk 
sound velocity v is a constant here, and the scattering is 
elastic, we can equivalently write Eq. as 



N 



T P h( £ )= J2 



fcn'e I, I 2 
~Z f-"n.n.' n 



(37) 



n.n'—0 



with the k nc given by Eq. i|ll|) . 



G. Landauer formula 



The charge current / and linear conductance 

G = lim — 

v-o V 



(38) 



are related to the electron transmission probability T e i(e) 
through the Landauer formula&i 



/>oo 

I= ^hj ^T el (e)[<(e)-<(e)]. 



(39) 



FIG. 5: Individual electron transmission probabilities |i n „/| 2 
as a function of energy for a circular right-angle bend with 
R/b = 1.20. Here A cl = % 2 h 2 /2mb 2 . 



Here rip(e) and rip 1 (e) are the Fermi distribution func- 
tions in the left (1) and right (r) leads, with chemical po- 
tentials /ii and /i r differing in proportion to the applied 
voltage V = (m — Mr)/e. 

Similarly, the thermal current I t h and conductance 



G 



ih 



lun — — 

AT^O AT 



(40) 



are determined by the phonon transmission probability 
Tph(e) according t oV ' 11 



1 f°° 

/th = 2^h I Tph(e) ^ (e) _ (*)] ■ 



(41) 



rig 1 (e) and rig r (e) are Bose distribution functions in the 
left and right leads, with temperatures T\ and T r . 



III. ELECTRON TRANSPORT 

We turn now to an application of our method to coher- 
ent electron transport through a circular right-angle bend 
with outer radius R and width 6, as shown schematically 
in Fig. In this case the curvature profile is 



k(X) 



X), 



(42) 



where O is the unit step function. The origin of the X 
coordinate is taken to be one of the locations where the 
straight and curved sections of the wire meet. 

As explained above, the numerical method we 
use to calculate t nn i (e) requires the matrix elements 
((pcmelVl^a'n'e') of the effective potential ij2"T|) in the un- 
perturbed scattering states i|10|) . However, matrix ele- 
ments of the second term in Eq. 1)21(1. which contain the 
curvature gradient 



1 



k>{X) = -[6(X)-8(X-¥)] 



(43) 



2 4 6 



el A 



10 12 14 



FIG. 6: Total electron transmission probability T e i. System 
parameters are the same as in Fig. [K] 



involves integrals of a delta function S times functions 
F(Q) of 0. Integrals of this form, involving products 
of generalized functions, have to be evaluated carefully, 
as we show in Appendix ^ Apart from this technical- 
ity, the application of our method to this geometry is 
straightforward. 

Mesoscopic charge transport through bent wires has 
already been studied extensively, 25 and we will only con- 
sider one case of this, namely R = 1.2 b. First we cal- 
culate the individual electron transmission probabilities 
|£nn'(e)| 2 from Eq. 1|33|1 for the lowest few channels n 
and n' . The results are given in Fig. [S] and are in excel- 
lent agreement with the mode-matching results of Sols 
and Macucci^i for the same value of R/b (see Fig. 2a 
of Ref. l3lh. Our result for |£n(e)| 2 also agrees qualita- 
tively with that calculated by Lin and Jaffe^ for a right- 
angle bend with a slightly larger curvature (see Fig. 8 of 
Ref.El. 

The total electron transmission probability T c i(e), de- 
fined in Eq. 1)35(1 . is presented in Fig. for the same 
curved wire. At energies given by 



n 2 A r 



with 



1,2,3, 



(44) 



where A e i = TT 2 h 2 /2mb 2 , additional channels in the wire 
become propagating and contribute to the transmission 
probability. The threshold energies (|44|) follows from 
Eq. i|llf> and Fig. EI The principal effect of the curvature 
in the wire is to soften the transitions at these thresholds. 



IV. PHONON TRANSPORT 

We turn now to the main emphasis of our work, the cal- 
culation of transmission probabilities for two-dimensional 
scalar phonons with energy e = hut and (bulk) sound ve- 
locity v to propagate through curved wires. We are not 
aware of any previous work on this problem. 

As before, we consider a circular right-angle bend with 
outer radius R and width &, as shown schematically in 
Fig. with curvature profile given by Eq. I|42|) . The 




FIG. 7: Individual phonon transmission probabilities \t nn i \ 
as a function of energy for a circular right-angle bend with 
R = 2b. Here A ph = nhv/b. 




FIG. 8: Individual phonon transmission probabilities |t n „/| 2 
as a function of energy for a circular right-angle bend with 
R = 1.001 b. 



method of solution is the same as that outlined in Sec. If HI 
and Appendix lAl except that the transverse parts of the 
unperturbed scattering states, defined in Eq. (|12f) . are 
now different. 

In Fig. we give the individual phonon transmission 
probabilities \t nn >(£)\ 2 , defined in Eq. (|33|) . for the lowest 
channels in a smoothly bent wire with R = 2 b. In Fig. [5] 
we do the same for a more tightly bent wire, with R = 
1.001 b (inner radius is 10~ 3 b). At energies given by 



e = ?iA c i with n = 1, 2, 3, 



(45) 



where A p h = nhv/b, additional channels in the wire be- 
come propagating and contribute to thermal transport. 

In both examples, transmission is nearly perfect in the 
low-energy e < A p h limit, where there is only a single 
propagating channel. At higher energies, scattering does 
occur. However it is mostly in the forward direction, 
and the fraction of transmitted energy T p h(e), defined in 
Eq. I|36|) , is essentially unchanged from that of a straight 
wire. 



FIG. 9: Dimensionless thermal conductance Gth/Gq as a 
function of temperature, for a 100 nm curved Si-like quantum 
wire, with outer radius R = 1.0016. Here G q = nk%T/6h, 
which is itself linearly proportional to temperature. 



ph 

FIG. 11: Individual phonon transmission probabilities \t nn /\ 2 
as a function of energy for a circular right-angle bend with 
R= 1.00001 b. 




0.99- 



'0 0.1 0.2 0.3 0.4 

T (K) 

FIG. 10: The solid curve is the same as Fig. H Dashed 
curve is the dimensionless thermal conductance for a straight 
Si-like wire with b = 100 nm. Thermal transport is hardly 
suppressed by the bending. 

In Fig. [3] we plot the thermal conductance Gth hi units 
of the "quantum" of thermal conductance 

G q ^^«0.95TpWK- 2 , (46) 
on 

for a curved wire of width 6 = 100 nm and outer radius 
R = 1.0016. The scalar phonon velocity is taken to be 
v = 8.5 x 10 5 cm s _1 , the longitudinal sound speed in Si. 
The thermal transport is hardly affected by the curvature 
in the wire, as can be seen in Fig. 1101 which compares 
an expanded plot of G t h/G q to that for a straight wire. 
The greatest suppression occurs near 200 mK and is only 
about 0.5% of the thermal conductance quantum. 

It is physically unrealistic to consider a 100 nm wire 
bent more sharply than R = 1.0016, because the in- 
ner radius of 0.10 nm in this case is already approach- 
ing atomic dimensions. However, for a wire of width 
6=10 (im and the same inner radius of curvature, we 
have R = 1.000016, the transmission probabilities for 



which are shown in Fig. 1111 The transmission prob- 
abilities when R = 1.000016 are similar to that for 
R = 1.001 6, shown previously in Fig. [HI as is the thermal 
conductance. We find in that in a 10 /^m Si-like quantum 
wire with R — 1.00001 6, the greatest suppression in G t h 
occurs near 2mK and is again only about 0.5% of the 
conductance quantum in magnitude. 



V. DISCUSSION 

We have introduced a general method to calculate the 
transmission of scalar waves appropriate for mesoscopic 
electron and phonon transport through a curved wire or 
waveguide. Applications to electron transport accurately 
reproduce results obtained by other methods. Phonon 
transport through curved wires is considered here for the 
first time. 

Our results demonstrate that curvature hardly sup- 
presses thermal transport, even for sharply bent wires, 
at least within the two-dimensional scalar phonon model 
considered. This behavior can, to some extent, be un- 
derstood by considering transport in the extreme long- 
and short-wavelength limits. In the long- wavelength, low 
energy limit, T p h — > 1, a consequence of the rigid-body 
nature of the wire. T p h also approaches unity for short 
wavelengths, because in this limit the phonons cannot 
sense the curvature. 

Because the phonon reflection probability always van- 
ishes in the long- wavelength limit, a simple perturbative 
(Born approximation) treatment is posssible at low en- 
ergies. For example, the energy-dependent n = trans- 
mission probability is 

|t 00 | 2 = 1 - kool 2 , (47) 
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where, to leading nontrivial order, 

roQ = fd 2 Rd 2 R' e^ x / h " 

2eb J 

x V a G (R,R',e)V R >e ieX '/ hv . (48) 

V R and V R i act on R and R', respectively. This result 
allows the low-temperature thermal transport though a 
variety of wire shapes to be addressed quite simply. Al- 
though an analogous perturbative expression can be de- 
rived for the electronic transmission probability as well, 
the form of the transverse part \n of the unperturbed 
scattering states, as dictated by the hard-wall bound- 
ary conditions, then leads to a divergence in the Born 
series^ consistent with the fact that T c i — * in the 
long- wavelength limit. 

We conclude with a brief discussion of the experimen- 
tal implications of our mesoscopic thermal transport re- 
sults, the charge-trasport case having already been dis- 
cussed in the literaturei^M Thermal transport in car- 
bon nanotubes has been studied experimentally by sev- 
eral groups and nanotubes would be interesting 
systems to use to investigate the effects of bending on 
transport. To apply our method of analysis to this sys- 
tem would require the consideration of scattering of clas- 
tic waves in a curved, hollow tube. Although they were 
obtained for scalar phonons in two-dimensional strips, 
our results do suggest that the effects of curvature will 
be small, if not completely negligible, in these systems. 
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APPENDIX A: INTEGRALS INVOLVING 
PRODUCTS OF DIRAC DELTA FUNCTIONS 
AND STEP FUNCTIONS 

Here we discuss integrals of the form 

/oo 
dxF[Q(x)]S(x), (Al) 
-oc 

where F is twice continuously diffcrcntiable, and 9 (a;) 
and S(x) are the unit step and Dirac delta functions, 



respectively. Integrals of the form ijAlfl . which involve 
products of generalized functions, depend sensitively on 
how 9(x) and S(x) are defined. It will be necessary to 
define 9(.t) and S(x) according to their appearance in 
this work. 

The step function Q(x) appearing in Eq. IjAlf) origi- 
nated from the curvature profile of Eq. (|42[1 used to de- 
scribe a circular segment of wire connected to a straight 
lead, and the delta function comes from its derivative 
with respect to arclength in Eq. (|43[) . which is required 
by the effective potential V of Eq. 1|21[) . Therefore we 
require Q(x) to be smooth (on some microscopic scale) 
and continuous, and S(x) to be related to it by 

^e(s) =<*(*). (A2) 

We also require, of course, that 

lim Q(x) = and lim Q(x) = 1. (A3) 

x — y — oo x — >oc 

The precise shape of O(x) near x = is immaterial, but 
with no loss of generality we can require that 9(0) = |. 

Integrals of the form IjAlJ) are now well defined. For 
example, 

/oo 
dxQ{x)5{x) = \, (A4) 
-oo 

as expected, but 

/oo 
dx[Q{x)] 2 5{x) = ±, (A5) 
-oo 

instead of j. These results are obtained by integrating 
by parts and using the behavior of 9(x) as x — > ±oo, not 
by evaluating 9(0) and [9(0)] 2 . More generally, 




dx [<d(x)] n 5{x) = -, (for n > 0) (A6) 



which is different from the naive value of [9(0)]" = (±) n , 
unless n = I. 
The reason why 

/oo 
dx F[e(x)] 5{x) ^ F[6(0)] (A7) 
-oo 

in some of these examples is because the delta function 
is distributed over a small but finite region of x, whereas 
9(x) and F[Q(x)\ are generally not slowly varying over 
that length scale. We conclude, therefore, that integrals 
of the form lAlf) appearing in the evaluation of matrix 
elements of V, have to be evaluated using intcgration-by- 
parts (or with an equivalent method). 
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